DYNAMIC MODULI AND LOCALIZED DAMAGE IN COMPOSITES 


L.S. Fu 

The Ohio State University 
Columbus, Ohio 43210 


The aim In nondestructive evaluation (NDE) by sound and ultrasound Is to 
extract material and fracture properties from pulses sent by and received at trans- 
ducers. The analysis Is usually given either In the time domain or the frequency 
domain. The Introduction of the size factor, the wavelength, In the model mechanics 
problems In NDE plays the Important role of relating the far-fleld measurable quan- 
tities such as phase velocities and attenuation to the near-field physical situation 
such as local geometric dimensions and elastic properties. This aspect of the anal- 
ysis thus allows “nondestructive" testing methods to be employed for "experimental" 
validation of the predictions by the theories (refs. 1 to 3). 

Solution to the elastic wave scattering due to a single embedded Inhomogeneity 
Is available by different methods that are appropriate at different frequency 
ranges. The methods that offer a solution In an analytic form and are useful for 
Inhomogeneous media with multiple components are the long-wave approximation, the 
polarization approach, and the extended method of equivalent Inclusion (refs. 4 
and 5). 

Several averaging schemes or theorems exist In the literature for finding the 
dynamic effective moduli and mass density. Some efforts concentrate on the average 
stress & and strain £ fields or on the average displacement field u,; others use 
a variational approach. These theories are appropriate mostly at RayleTgh or long- 
wave limits and do not exhibit dispersive effects. Dispersiveness and attenuation 
are Important In evaluating dynamic material properties. 

This study dealt first with the scatter of elastic waves due to a thin, flat 
ellipsoidal Inhomogeneity, either penny shaped or elliptical. An average theorem 
appropriate for dynamic effective mass density and effective moduli was developed 
via a self-consistent scheme. Effective material properties of two- component media 
consisting of randomly distributed spheres are given here as a special case. 


PRELIMINARIES 

Displacement Field due to Presence of Mismatch 

The Inhomogeneous media considered In this paper are assumed to consist of a 
homogeneous matrix of elastic moduli £? and mass density p° and a distribu- 
tion of inhomogeneities with moduli C^ 1 ") and mass density p( r ) occupying 
regions ft r , r = 1,2,..., n (fig. 1). The total displacement field u can be 
separated Into two parts. 


u = uO) + u( m ) ( 1 ) 

where the superscripts (i) and (m) denote "Incident" and "mismatch," respectively. 

It Is clear that when no mismatch components (i.e., inhomogeneities) are present. 
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the total displacement field Is entirely the same as the Incident displacement wave 
field. On the other hand, If there Is no Incident wave field, the only field that 
exists Is the null field. 


Elgenstralns 


The elgenstralns e,*, which are also termed transformation strains, are 
defined as the part of the total strain & that must be subtracted before the 
remaining part can be related to stresses" a through Hooke's law 


C = £ 

rs rs 


- c 


rs 


( 2 ) 


e rs = C jkrs°jk (3) 

where C 7 1 are the elastic compliances. The method of equivalent Inclusion Is a 
method That allows the Inhomogeneous media to be replaced by media of homogeneous 
matrix effective moduli with distributed transformation strains In the regions 
originally occupied by Inhomogeneities; hence 

0 In Q 

* ° (4) 

e rs in G r , r = 1 , 2 ,.., n 

For the two problems to be equivalent, the transformation strains must give rise to 
a field that Is exactly the same as the "mismatch 11 field ,u( m ). This leads to 
equivalence conditions that ensure that identical field quantities at any given 
point will be obtained In the two problems* (ref. 4): 

4 C jkr S u r?s< f) '♦ c jkrsv! 1) < ? > ' ' 4 C jkrs u r!s< f > 1nft < 5 > 

A / o6> 2 u^ m) (r) + Cj krs e*^(r) = -Apu 2 Uj^(r) In Q (6) 



Two types of transformation strains, or elgenstralns, arise In elastodynamlc 
situations due to the mismatch In elastic moduli A ,c and mass density Ap. It Is 
often convenient and useful to define associated quantities such as 


*(11 

n jk = C jkrs e rs 


(7) 


" C jkrs e rs,k 

where mj k and irj are moment density tensor and equivalent force or 
elgenforce, respectively. 


(8) 


+The conditions (eqs. (5) and ( 6 )) are similar to those of Willis (1980) 
and those of Mura, Proc. Int. Conf. on Mechanical Behavior of Materials . 5, Society 
of Materials Science, Japan, 12-18 (1972). 
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Volume Average and Time Average 


The volume average and time average of a field quantity, say F(r,t), are 
denoted by using brackets <> v and <> T , respectively, and are defined as 


<F(l,t)> = l/F(v,t)dV 

(9) 

<F(l,t)> = j J*F( r ,t)dT 

(10) 


where V and T are volume and time period, respectively, and v Is velocity 
field. 


Volume Integrals of Ellipsoid Associated with Inhomogeneous 

Helmholtz Equation 

Volume integrals of an ellipsoid associated with the Integration of the 
Inhomogeneous Helmholtz equation are used In this work. The Inhomogeneous scalar 
Helmholtz equation takes the form 

V 2 $ + k 2 $ = -4iry(r) (11) 

where y(l) is the source distribution or density function and V 2 and k 
are the Laplacian and wave number, respectively. A particular solution to 
equation (11) is 

$(r) = £ , exp(ikR)dV', R = |r - r 1 | (12) 

in which (4irR)-! exp(lkR) is the steady-state scalar wave Green's function and 
Q Is the region where the source Is distributed. The source distribution function 
Y(r) can be expanded in basic functions or polynomial form, depending on the geom- 
etry of the volumetric region «. For an ellipsoidal region, the choice of using 
a polynomial expansion separates this work from other theories of elastic wave 
scattering : 


Y(r') = (x')Hr)' i (z') v (13) 

In which p, and v are integers and k is either longitudinal or transverse 
wave number. 


For elastic wave scattering in an Isotropic elastic matrix, two types of 
volume Integrals and their derivatives must be evaluated: The ^-Integral s are 

given by 


<k(r) = R _1 exp(1aR)dV' 

(14a) 

*k ( r) = J £ x'R -1 exp(1aR)dV 

(14b) 

s (l) = J £ x£x^. . .x^R -1 exp( iaR)dV 1 

(14c) 
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(15a) 


♦ . p (r) = *(£.) 

v P (I) ■ hr *k<t> < 15b > 


2 2 2 2 

where a = u> /V. = pw /(\ + 2y) . The other type, the ^-Integrals, are 
obtained by replacing a with 6 In equations (12) and (13), where 
6 2 = pw 2 /(y). Details of the Integration are given In reference 6. In this 
work, only the limiting values at r -» 0 and r -» <» are of Interest. 


ISOLATED, FLAT ELLIPSOIDAL INCLUSION 
Formulation 

Consider the physical problem of an Isolated Inhomogeneity embedded In an 
Infinite elastic solid that Is subjected to a plane time-harmonic Incident wave 
field as depicted In figure 1. Replacing the Inhomogeneity with the same material 
as that of the surrounding medium, with moduli Cj|< rs and mass density p, and 
Including In this region a distribution of elgenstralns and elgenforces, the 
physical problem Is now replaced by the equivalent Inclusion problem. 

The total field Is now obtained as the superposition of the Incident field and 
the field Induced by the presence of the mismatches In moduli and In mass density 

*( 1 ) * 
written In terms of elgenstralns and elgenforces ir^ 

F = fO) + F( m ) (16) 

where £ denotes either the displacement field uj, the strain field e^j, or 
the stress field 

For uniform distributions of elgenstralns and elgenforces, the fields can be 
obtained as 


u i m> - C Jkrs4s 1)s 3m ,k' <£,) 

(17) 

e (m) = [ (m) (m)l / 

mn [ u m,n + u n,mJ/ 1 

(18) 

»<"■> . C c (m) 

pq pqmn mn 

(19) 


where a comma denotes partial differentiation and 

V £) '{v i - I,,dv ' (20) 

In which gj m Is the spatial part of the steady-state elastic wave Green's func- 
tion and Q. Is the region occupied by the Inhomogeneity. The Integrals Sj m 
and their derivatives must be evaluated for the regions r > r 1 and r < r 1 
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(ref. 6 ). The solution form represented In equations (14) to (18) gives the fields 
Inside and outside an Isolated Inhomogeneity of arbitrary shape. 


Far~F1eld Scattered Quantities 

Let the Incident displacement field be longitudinal and of frequency w and 
amplitude u 0 : 


.0) 


= Vj 


exp( 1 ax^k^ - Iwt) 


(21) 


where I 2 = -1 and q^ Is the unit vector In the normal direction of the plane 
time-harmonic Incident wave and Is wave vector. For a linear Isotropic 
medium, the spatial part of the free-space Green's function Is well known. Substi- 
tuting gj m (r - r') In equations (18) and ( 20 ) and using the limiting concept 

11 m a 3 irj = Aj constants ( 22 ) 

a 3^0 

*/ 1 ) 

11m a 3 e ij = B 1j constants (23) 

v° 


the scattered displacement u l ( 1 s )(r,t) from a thin, elliptical flat crack can 
easily be obtained as 

m - m — ' 


(aa,) 3 u (aa n ) 3 u 

v 1 ' 0 v 1 ' 0 


r-*» 


where 


H m = (a 


m 


= j^(CG m exp 1ar)y/ar + (DH m exp 1 J3 r (3 r exp(-lcot) 

G m - -<v*i> [-v/j * <’ - * 2 <“ 2 / 62 >Vk l 3 B ki] 

/a,) [-<BA.) 3 ( Vj - 5 mj» A j - 2 < 8 / “> 2 l k B km * 2 < B/ “>VkVkJ 

C = [j 0 *“ r e' * ■M“ r e ) ] / / 3 ° s1ri “ r e / <“ r e> 3 " cos “ r e / *“ r e )2 
D = * J 2 (B r e )j 1 3 - s1n Br/tBr^ 3 - cos Br e /(Br e ) 2 


? 7 ? 

r e = Vi’ a 3 = 0 


(24) 


(25) 


(26) 


In which m,j,k = 1,2,3 and (a-|,a 2 )» S-m, a * and B denote the semiaxes of the 
flat ellipsoid, direction cosines of scattered displacements, longitudinal wave 
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number, and shear wave number, respectively. Also, Aj and Bj|< are the reduced 
nondimensional forms of Aj and Bj|<, respectively, defined as follows: 


Aj = Aj/( A/ou 2 U 0 ) , A p = p' - p 

Bjk = -B jk/(^ a ^o) 


(27) 

(28) 


Expressions for the differential cross section dP(<o)/dQ and the total 
cross section P(w) can be obtained as (ref. 5): 


= <j L (e,<t>) + (a/B)o T (e,<t») 


P( 


<j>) = ( © , <J>) + (a/B)a (0,<j>)j 


dn 


(29) 

(30) 


where dn Is the differential element of solid angle and 


2 L 


aV(e,«t,) = (aa 1 )°[CG m ][CG m ] 


(31) 


B 2 d T (e,4>) = (aa 1 ) 6 [DH m ][DH m ] (32) 

In which the overbars denote complex conjugates. The constants A? and Bj|< 
must be evaluated from the equivalence conditions (eqs. (5) and (6)) with the use 
of the limiting concepts In equations (21) and (22) and of the Integration method 
developed In reference 6. 


Determination of Aj and Bj|< 

In equations (21) to (25) the scattered displacement field Is given In terms 
of the "reduced 1 forms of the elgenforces and elgenstralns (l.e., Aj and Bj|<). 
These constants must In turn be determined from the equivalence conditions. By 
writing the Incident wave field In a Taylor series, the governing simultaneous 
algebraic equations can be easily obtained. Since fjj[0] and F m jj[0] vanish 
automatically, these governing equations become uncoupled and lead to a three-by- 
three system for Aj and a six-by-six system for Bj|<. For a linear elastic 
medium, they are 

Apa> 2 u 0 f j $ [ 0 ]A j + As = -q s (33) 

{A\6 s tD mm j j [0] + 2AvtD s -tj kfO] IBjk + (* & st B mm + 2yB s t) 

= — ( AX& s -tQni^nn + 2Ayq s q-t) (34) 

where the subscripts s,t,m,j,k = 1,2,3 denote sum from 1 to 3 and 

4irpo)2f j S ( r) = -B 2 4>5j S + <k*mj ~ 4>.mj ( 35 ) 
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(36) 


4irp« 2 D s tjk(l) = 2y(t|», s tjk + ‘tMtjk) ~ v3 2 (4>» jt 5 ks + 4>,js 6 kt) - Xa2 t,mn 4 jk 
In which 

AX = -X, Ay = -y, A p = -p 

The $ - and ^-Integrals and their derivatives are evaluated by the method 
suggested In reference 6. Retaining terms up to aa-j or Ba-j of the fourth 
order, the constants are obtained as 



Aj = -qj/{u 0 w 2 Apf j [0] + 1), no sum on j 

(37) 


fjCO] = (fn[0], f 22 [0], f 33 [0]) 

(38) 


{B*} = [bij]-l [Cj],1,j = 1,2,3 

(39) 

(Blj> 

= -qiqj/[l + 5(<h[0] + 4> j [ 0 ] ) J , no sum on 1,j 



1 * j; 1,J = 1,2,3; 5 = 1/4* 

(40) 

where 

In equation (37) B* = B*-|, Bj? « B 22 » ^3 = ®33» and 



C ^ = (AX + 2 Ay) ^1 + q 2 ^, j = 1,2,3 

(41a) 


t»H = (X + 2y) + Xfi|»,jj[0] + 2 ( X + 2y)5<t>,-|-| [0] + 2yCi|», -j -| [0] 

(41b) 


b 12 = X + XO,jj[0] + 2yf»i», 11 [0] + 25gX<(>, 22 C 0 ) 

( 41 c ) 


b 13 = x + + 2yf»|»,-| -j [0] + 2X54>, 33 [0] 

( 41 d) 


b 21 = X + XC«|r,jj[0] + 2'|iC4 r , 22 C 0 ] + 2^X<|>, ^ [0] 

( 41 e) 


b 22 = (X + 2y) + Xfi|/,jj[0] + 2yfi|f, 22 [0] + 2?(X + 2y )4>, 22 [ 0 ] 

( 41 f ) 


b 23 = X + Xf*,jj[0] + 2y£V, 22 [0] + 25X4., 33 [0] 

( 41 g) 


b 31 = X + XCf,jj[0] + 2y C'P » 33 C 0 ] + 2X54>, 11 [0] 

( 41 h) 


b 32 = x * K'k.jjtO] + 2yC«j;, 33 [0] + 2X^4>, 22 [0] 

(411) 


b 33 = ( x + 2 ^) + ^,jj[0] + 2yCi|r, 33 [0] + 25 <t>, 33 (X + 2y)[0] 

(41 j) 

Note 
<t>3 * 

that b^j * b<^. In equation (38), <t>-| = 4>-)-][0], <t> 2 = 4>22t°]* and 
<t> 33 [0] . The f- and <t>-f unctions are given as 



4irp« 2 f j s [0 ] = -B 2 4>[0]6j s + xp, j $ [O] - 4> , j s [ 0 ] 

(42) 


B 2 4>[0] = 1 ra 1 a 2 B 2 Jl 0 - [(3a-, ) 2 /l&] ^ + 1 (8/3)3> 

(43) 


v / 

■- [(* a ? a 2 04 )/ 12 ] 'l 

(44) 
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In which 


<l>, 22 [0] = "[( ira l a 2 134 )/ 12 ] ! 1 


I o = 






= (* , » d * — = b - 4"l F( 0* k) 

J 0 (a; + *jA(«) a l L k k J 

T _ f jLj jH 2_ [~F(e, 

2 ’ J n fa? + *W> 3l L k2 


F = 


| o - .".Ar 1 • 1 (' ■ 


kl Ej eAl 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


and as a g -> 0, e * ir/2, k 2 -► (l - a 2 /a 2 ) and k' 2 = (1 -.k 2 ) -» a 2 /a 2 , If 

a-| > a 2 . If a 3 ■* 0 and a-| = a 2 , we have I 0 = ir/a-j and I-, = I 2 = ir/2a-| . 
The ijr-f unctions are obtained by replacing 6 with a In the <t>-f unctions 

(fig. 2 ). 


Numerical Calculations and Graphical Displays 

It Is clear from equations (37) to (40) that the uniformly distributed eigen- 
strains and elgenforces In their reduced forms, B-m and Aj , respectively, depend 
only on the characteristics of the Incident wave field and the geometric factors of 
the Inhomogeneity for a given material system. It Is observed from equations (37) 
to (41) that and B^ as functions of dimensionless wave number aa-| or aa 2 

would exhibit large peak values at certain Incident wave frequencies when b^ -» 0. 
The values of these critical frequencies depend only on the matrix elastic moduli, 
the crack dimensions, and the measurement direction. For a given aspect ratio a 2 /a 
and a measurement direction, the difference In frequencies at subsequent peak values 
Is proportional to a-] , the largest dimension of the Inclusion. 

Computational data of elastic wave scattering due to a flat embedded Inhomoge- 
neity In any given Isotropic material system can be obtained by employing equations 
(24) to (26) and (33) to (50). Scattered displacement amplitudes can easily be 
obtained for any given aspect ratio. In figures 3 and 4, computational data for the 
back-scattered situation are displayed for a tungsten disk In an aluminum matrix 
for aspect ratios a 2 /a-| = 0.01 and a 2 /a '\ = 0.50. In figures 5 to 8, critical 
frequencies were observed In bU, but no critical frequencies were Identified for 
A? (j = 1,2,3). The critical frequencies for C|G m | and D|H m | are, however, 
clearly Identified. The position along the aa-| axis at which the first critical 
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frequency occurred depended on the aspect ratio and matrix elastic moduli. A suf- 
ficiently small Increment In aa] must be used In order not to miss any peak 
values. Since the solution form given In equations (24) to (26) and (33) to (50) 

Is analytic In frequency, this can easily be achieved. Most of the scattered 
energy Is carried by the transverse components of the scattered displacement (l.e., 
I DH m | > I CG m | , figs. 3 and 4). 


DYNAMIC MODULI AND DAMAGE OF COMPOSITES 

Consider the problem of an Inhomogeneous medium as Illustrated In figure 9, 
under a plane time-harmonic Incident wave field. The true composite thus occupies 
the whole region and possesses effective moduli JT* and mass density p*. To 
determine the effective moduli and mass density, an average strain energy and 
kinetic energy were used. The effective properties were found to depend on a 
fourth-rank tensor A, and a second-rank tensor J). A self-consistent scheme was 
then developed for determining these tensors. 


Average Theorem 

To determine the effective moduli and mass density, the following definitions 
were used: 


<S> = £*<£> 

(51) 

<et> - £*<o<c> 

(52) 

<pv> = p*<v> 

(53) 

<pV‘V> = p*< V><V> 

(54) 


where a, &, v are the stress, strain, and velocity fields, the ~ denotes a 
tensorlal quantity, and the angular brackets O denote the volume average of a 
field quantity (eq. (9)). The left and the right sides of equations (52) and (54) 
can be shown to be equivalent under the so-called Hill's condition (ref. 7), <£><v> 

= <de>. From equations (53) and (54), It Is clear that the kinetic energy per unit 
volume of the effective medium can be made equal to that of the physical medium if, 
and only If, a frequency-dependent mass density is defined. This Is the same as 
requiring, by again using the Hill's conditions <p,v«v)> = <pv><()C>» the avera 9 e 
linear momentum per unit volume to be the same as the effective linear momentum per 
unit volume (eq. (53)). These conditions (eqs. (53) and (54)) are not met If the 
static definition of effective mass density Is used. 

Let f r denote the volume fraction of the r^h inclusion material. Then 
the volume averages of the stress and velocity fields a and u, respectively, are 

<s> - E f r e <r> <55) 

r=0 

<u> = (56) 

r=0 r 
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( 57 ) 


where a( r ) = (r) = 0,1,2, ...,n and for a time-harmonic situation 

u = v = -iwu 

Equations (52) and (54) can be rearranged, by using equations (55) and 
(56), as 


<s£>- c °<sXs> * E f r £, (r) e, (r) < £ > 


(58) 


r=l 


</>v-v> = p° <v><v> + £ f r Ap (r) v (r) <v> 

r=l r 


(59) 


Consider now the case where the solution form possesses a linear relation between 
the velocity and strain fields In the r th component and the average velocity and 
strain fields of the effective medium, that is, 

,< r ),A< r ><c> (60) 

y(r) „D( r ><v> (61) 

where A( r ) and JD( r ) are tensors of fourth and second rank that must be 
determined with a suitable scheme. Substituting equations (60) and (61) In 
equations (58) and (59) leads to 


<£&> ■ 

< P vv> * 


C° + £ f AC (r) A (r) 


r=l 
n 


■*If ^ (r) 0 <r) 


r=l 


<^Xe> 

<y><v> 


Comparing these equations with equations (52) and (54) gives the following 
expressions for the dynamic effective moduli and mass density: 


C* . c° ♦ £; f AC (r) A (r) 
r=l r 


»• = p" • Z f, ^ <r) 0 (r) 

r=l r 


D (r > » 0<;>/3 


(62) 

(63) 


(64) 


( 65 ) 


Note that the tensors /\( r ) and o( r ) are frequency dependent and replace the 
static expressions when frequency approaches zero and when proper care Is taken. 
The assumption of a general linear dependence between v( r ) and <v> (eq. (61) 
must be specialized such that the second-rank tensor D, will degenerate Into a 
scalar. This specialization Is automatic for randomly distributed spheres where 

D mj * D V)« wh1ch 0 * D ll = d 22 = d 33- 
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To determine the explicit form of the tensors A( p ) and JD( r ), the strain 
and velocity fields In the r^h component are determined by using the method of 
equivalent Inclusion as presented In the previous section. 


Self-Consistent Scheme for Determining Effective Properties 

The tensor fields A( r ) and D^ r ) are of ranks four and two, respectively, 
and are functions of wav£ numbers, geometric properties, and effective and Inclusion 
material properties. Let the average strain and velocity be the same as those 
derived from the Incident wave field. Then the governing conditions for determining 
A( r ) and D,( r ) for the r* h Inclusion are simply obtained by rewriting the 
equl valence~cond 1 t 1 ons for an effective medium with mass density p* and moduli 
C*: 

C(r)r c O) + e (m )1 = C *r e n) + e (m) - e *l In Q r ( 66 ) 

p( r ) w 2|yi) + u,(m)| = p*w 2 £uO) + u( m ) - ttJ In G r (67) 

In which the superscripts (m) and ( 1 ) denote mismatch and Incident, respec- 
tively. Slightly different approaches for finding effective moduli for heterogene- 
ous materials that apply to static cases with different constituents and situations 
are given in references 7 and 8 . The strengths and dynamic responses of composites 
have been Investigated (refs. 9 and 10) and reviewed (refs. 11 and 12). 


Effective Properties of Two-Component Media: Randomly Distributed Spheres 

Let the moduli, mass density for the matr1x f Inclusion material particles, and 
the effective medium by denoted by Cj|< rs , p°; Cj|< rs , p 1 ; and Cj|< rs and p * > 
respectively. By using the elastodynamlc solution for a single ellipsoidal Inhomo- 
geneity, the displacement and strain fields Inside an Inhomogeneity are found, when 
the average displacements are made equal to the plane time-harmonic Incident wave 
field: 



= D(r)u^) 


In S2 r 

( 68 ) 

e,( m >(r) 

= A(r) G n) 

— * r^y 

- £(l)y, (1) 

In « r 

(69) 


Employing the volume-averaging process as described by equation (9) and substituting 
In equations (58) and (59), the effective properties are easily defined as follows: 

P* = p + f ApD, 0 = 3 D jj (70) 

C* = C + f ACA (71) 

where f Is the volume fraction of Inclusion material. The tensor fields £ 
and A are ~ 

D mj = '< f mj (l) ^/i f MJ t0 ] * 4ir(p ' “ no sum on M ’ J ( 72 ) 
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A mnpq = ^mjk.n^ + F njk,n/-^ S jkpq/ 2p “ 


(73) 


In which the tensors f and F are defined as 


mj' 


mj 


* ♦•mj - 

Jr 74) 

•1 s mj “ 2 "♦•ml * 2 "*’mlj 

(75) 


The and, by Implication, the ^-Integrals given In equations (14) and 
(15), etc., are the volume Integrals associated with the Inhomogeneous Helmholtz 
equation. They can be carried out for an ellipsoidal region by expanding (exp 
1kR)/R In Taylor series expansions with respect to r 1 for r > r' and with 
respect to r for r < r*. Here k can be either a or 0. Details are given 
In reference 6. This type of expansion for the Integrand Is particularly useful in 
determining the coefficients of a "polynomial” distribution of irj and 

The fourth-rank tensor S^ pq Is the connecting tensor between the 
elgenstralns and the applied strains, that Is, 


*(D_ S « (a) 

jk " jkpq pq 


jkpq 


= S jkpq (£ ’ 


<*>) 


(76) 


for the case of uniform elgenstralns and elgenforces. In developing these 
expressions, the volume average of the <j>-1ntegrals must be evaluated. Finally, 
note that P * and C,* are complex, where the real and Imaginary parts are 
associated with the Velocity and attenuation, respectively. 


Example: Spherical Inclusion Materials 


Let the spherical Inclusion materials of radius a be randomly distributed 
over the whole volume of the matrix. If the matrix and the Inhomogeneities are 
Isotropic, the effective medium Is also Isotropic. It Is easy to show that 


D 


mj = 4 mj 


-< f 33<a> 

f 33 [0] + 4ir(*>' - p*)u> 2 


6 ,D 
mj 


(77) 


and 


S jkpq = S kjpq = S jkqp = S pqjk 
S llll = S 2222 3 S 3333 = C 1 

S 2323 = S 1 31 3 3 S 1 21 2 = C 3 

S 1 1 22 = S 1 1 33 3 S 2233 3 C 2 


(78) 
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where 



C 3 = l 2F 122,l^ 0 ^ + " v *^| 

c 2 - GF m,i^ + < G + 1 ) F 122,lt°] + H 
c i* ■ F m ,i [0] + 2 GF i22,i^ - F 

F = -(X* + 2y*)/G 

G = (X 1 - X*)/[(X' - X*) + 2(y' - y*)] 

H = X*/G 

Following the theory developed In the previous sections, the effective moduli 
and mass density are found to be 


p* = P + f ApD (79) 

X* = X + f[AX(A ini + 2A n22 ) + 2 A V A 1122 ] (80) 

Vi* = y + f a V>(A 1212 + Ai 2 2i ) ( 81 ) 

K* = K + f[(A ini <■ 2A 1122 )AX + (2/3)(3A 1122 t A 1212 + A^^jAy] (82) 


Clearly, the velocities are dispersive. At frequencies above the Rayleigh limit, 
this phenomenon Is pronounced. From figures 10 to 12, the bulk moduli, shear 
moduli, and longitudinal velocities are shown as functions of the volume concentra- 
tion of spherical Inclusion materials for aluminum spheres In germanium for differ- 
ent dimensionless wave numbers aa. For a given fixed concentration, the moduli 
K* and y* and the velocities V|_ and vj are Increased as the dimensionless 
wave number aa Is Increased. The dispersiveness of effective shear modulus Is 
minimal and that of effective bulk modulus Is more pronounced (figs. 12 and 13). 

As an example of detecting localized damage by void nucleatlon, let all small 
voids be locally nucleated within a localized small region n of radius R 
(fig. 14). The effective moduli of this composite can therefore be obtained from 
equations (64) and (65). If void nucleatlon outside the region Q can be Ignored, 
the scattering of the composite sphere can easily be obtained. By using the com- 
puter program developed In reference 17 for that of a single spherical Inclusion, 
the scattering cross section for a composite sphere consisting of small voids In 
titanium Is displayed as a function of dimensionless wave number for different con- 
centrations of voids (fig. 15). As the volume fraction of voids Inside n Is 
changed, the effective properties p* , X*, and y* are also changed. Hence, the 
attenuation effect Is pronounced as the concentration of voids Is Increased. The 
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scattering cross section, which Is essentially proportional to the attenuation 
(ref. 1) Increases with Increasing concentration f r . It appears that these 
curves can be used to locate and calibrate porosity In a structural component. 


CONCLUDING REMARKS 

The velocity and attenuation of ultrasonic waves In two-phase media were 
studied by using a self-consistent averaging scheme. The effective medium had to 
possess the same strain and kinetic energy as the physical medium. The concept of 
volume averaging for physical quantities was employed. The solution depended on 
the scattering of a single Inhomogeneity. The theory Is general and can be applied 
to any multicomponent material system. Since the scattering of an ellipsoidal 
Inhomogeneity Is known, the average theorem presented In this report can be used to 
study the velocity and attenuation of distributed Inhomogeneities of shapes such as 
disks and short fibers. The orientation of these Inhomogeneities as well as their 
sizes, as In the spherical geometry, will necessarily Induce anisotropy In the 
effective medium. The scheme developed herein was not compared with test data. 

The use of a self-consistent scheme In determining static material properties for 
composites has been explored (ref. 13). 

Results for randomly distributed spherical Inclusions of radius a are pre- 
sented. Effective moduli and mass density were found to be dispersive. A simple 
model of localized damage was studied. It Is well known that porosity Is directly 
related to the strength of rocks and ceramics. Therefore the theoretical study of 
velocity and attenuation In two-phase media may be a viable means for data analysis 
In ultrasonic evaluation of dynamic material properties for composite bodies 
(ref. 14) and polycrystals (ref. 15). As In determining the static properties of 
materials, the problem may be more In the mechanics than In the manufacturing 
(ref. 16). Deformation processes such as rolling, sheet metal forming, and drawing 
often Involve large plastic flow, moisture absorption, and thermal cycling. Resid- 
ual stresses and anisotropy are Introduced Into the material and limit the amount 
of deformation to fracture with a directional dependence. Theoretical development 
In acoustoelasticity appears to be lacking (ref. 18). 
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Figure 3. - Ic^l versus aa-| for back-scattered d 
32/3! = 1; a! = 0.5; $ = 25°; 6 = 30°; 



Figure 4. - IDI^I versus aa-| for back-scattered c 
a 2 /a l 0.5; $ = 25°; 6 = 30°; % = -l m . 


2 




Figure 5. - |Af | versus aa] for back-scattered data for tungsten disk in aluminum matrix. 
32/3] = 1; a] = 1.0; <(> = 25°; 0 = 30°; Cfo = -t m . 



Figure 6. - |B]]| versus aa] for back-scattered data for tungsten disk in aluminum matrix. 
a 2 /a] = 1; a] = 1.0; <t> = 25°; 0 = 30°; = - 8 . m . 



Figure 7. - | B 13 1 versus aa] for back-scattered data for tungsten disk in aluminum matrix. 
33 / 3 ] = 1; a] = 1.0; 4* = 25°; © = 30°; = -!>m. 
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Figure 8. - Ig^ versus aa-| for back-scattered data for tungsten disk in aluminum matrix. 
32/3! = 1; 3! = 1.0; 4> = 25°; 9 = 30°; cfo = -l m . 
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Figure 11. - Longitudinal wave velocity versus concentration: aluminum spheres in germanium. 



DIMENSIONLESS WAVENUMBER aA 

Figure 12. - Effective shear modulus versus wave number: aluminum spheres in germanium. 
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DIMENSIONLESS WAVENUMBER aA 

Figure 13. - Effective bulk modulus versus wave number: aluminum spheres in germanium. 
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Figure 14. - Schematic diagram of localized damage in infinite solid. 
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